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Abstract 


In this paper an exact formulation of the plane elasticity problem 
for a hollow cylinder or a disk containing a radial crack is given. The 
crack may be an external edge crack, an internal edge crack, or an 
embedded crack. It is assuned that on the crack surfaces the shear 
traction is zero and the normal traction is an arbitrary function of r. 
For various crack geometries and radius ratios, the numerical results 
are obtained for a uniform crack surface pressure, for a uniform pressure 
acting on the inside wall of the cylinder, and for a rotating disk. 


1. Introduction 

In this paper we reconsider the plane problem of a hollow cylinder 
containing an arbitrarily oriented radial crack and subjected to arbitrary 
normal tractions on the crack surfaces. Aside from its direct applica- 
tions to plane problems such as cracked rings and rotating disks, the 
problem has potentially important applications to pressure vessels and 
piping containing a relatively long part-through crack in a meridional 
plane. In this latter group of three dimensional problems the plane 
strain solution provides an upper bound for the stress intensity factor 
in the mid-region of the longitudinal part-through crack in the cylinder 
wall. Because of the importance of the crack geometry, the problem has 
been studied rather extensively (see, for example, [1-7]). After 
Bueckner's early work in which the technique of complex potentials was 
used [1], most of the subsequent analytical studies were based on the 
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mapping technique originally developed by Bowie [2-4]. In [1] the 
problem of a rotating disk having a radial edge crack on its inner boun- 
dary is considered. The same crack geometry under uniform tensile 
tractions applied to the outer boundary of the ring is studied in [2]. 

The problem of multiple cracks located on the inner or outer boundary 
of the ring is considered in [3]. The problem of a curved beam or a 
ring segment containing a radial crack is studied in [4]. In [5] the 
problem of a hollow cylinder containing an external edge crack is con- 
sidered for various loading conditions where a numerical technique 
similar to finite difference approximations is used for the solution. 

In [6] the ring problem with an inner edge crack is solved by using the 
method of weight functions. Needless to say, because of its plane 
geometry the problem is also ideally suited for a finite element treat- 
ment (see, for example, [7]). 

There are, however, still sufficient reasons for reconsidering the 
problem. First, it may be observed that the existing solutions in one 
form or another are all restricted with regard to loading and/or geometry. 
Secondly, the published results are not, in most cases, sufficiently 
accurate in their entire domain of coverage. Third, it is desirable to 
have an exact analytical formulation and a reasonably simple and reliable 
method of solution for the problem which can be used by other investiga- 
tors to obtain results to any desired degree of accuracy for particular 
crack geometries and loading conditions. Finally, the embedded crack 
problem has not been considered in any of the previous studies. In addi- 
tion to its direct application to such problems as radial weld defects 
in longitudinal seams, the solution may be important in studying the 
crack problem in cylinders under self-equilibrating residual stresses 
(e.g., hollow tempered glass cylinders undergoing "static fatigue"). 

In this paper we consider the plane strain or the generalized plane 
stress problem for a concentric circular cylinder or ring containing a 
radial crack. It is assumed that the external loads are symmetric with 
respect to the plane of the crack and the basic ring problem without 
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the crack has been solved. Thus the main problem of interest is the 
Mode I crack problem in which the self-equilibrating arbitrary normal 
crack surface tractions are the only external loads (Figure la). 

It should be pointed out that in the special cases of solid disk 
(Figure la, a=0, b<») and the infinite plane with a circular hole 
(0<a<», b=®), the problem may readily be formulated in terms of a singu- 
lar integral equation by using the basic dislocation solutions given 
in [8] as the Green's functions. 


2. Formulation of the Problem 

The geometry and loading for the problem under consideration is 
shown in Figure la. In crack problems involving finite geometries, 
almost invariably the basic technique to formulate the problem is to 
express the solution as the sum of certain number of suitable solutions 
satisfying the differential equations of the problem in such a way that 
at least one of the solutions contains the main features of the general 
crack problem and the total solution contains sufficient number of arbit- 
rary functions or sets of arbitrary constants to account for all boun- 
dary conditions. In this problem the solution is expressed as the sum 
of a crack or dislocation solution for an infinite plane and the general 
solution for a circular ring. In order to facilitate the application 
of the boundary conditions, the dislocation solution is expressed in 
polar rather than in rectangular coordinates. Let, then, the stress 
state in the cracked ring problem be of the following form: 

j(r,e) = ct-, ( r , 0 ) + a 2ij .(r,e), (i , j = r,e) (1) 


where • and (i,j = r,e) respectively refer to an infinite plane 

containing an edge dislocation or a crack along 9=0 line and to a con- 
centric circular cylinder. 


a\ The infinite plane -solution. 

Consider an infinite plane with an edge dislocation having a Burger's 
vector by = -f located at the point r=t, 9=0 (Figure lb). The plane 
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elasticity problem may be solved by assuming that 

a lrQ = 0, 0 <_ r < ®, 0 = 0, 0 = it (2) 

[u l0 (r,+O) - u-] 0 (r,-O)] = fS(r-t), 0 < r < ». (3) 

Referring to [8] the Airy stress function of the problem may be expressed 
as (Figure lb) 


■h(r,e) = r, log r, cos 8, 

= “ ^ fcy f(r cose-t) log (r 2 +t 2 -2rt cose), 


(4) 


where u is the shear modulus, k=3-4v for plane strain and <= (3-v)/(l+v) 
for the generalized plane stress, v being the Poisson's ratio. From 
(4) the stresses are obtained to be 

-| 3i|>-j -| 3 2 i/)-| 

a lrr (r ’ 0) = r Tf“ + r 2 * ie 2 * 


2uf r rcose - t - 2t sin 2 e 
” tt(k+ 1 ) L r^+t 2 - 2rt cose ~ 

2t 2 sin 2 e(r cose - t) -i 
( r z+t^ - 2rt cose) z J ’ 

a l ^1 

a lre^ r ’ 9 ^ = “ 3r ^ r le“ ^ 

2yf r sine(2t cose-r) 

“ ir(ic+) ) L r^+t* - 2rt cose 


(5) 


2t sine(r cose-t)(r-t cose) 

(r 2 +t* - 2rt cose) 2 -* ’ 

3 2 4»-j 

a lee^ r ’ 8 ^ “ Tr 2 ” 

_ _ 2yf r 2cose(r-t cose) + r cose-t 
tt(k+ 1 ) L r z +t z - 2rt cose 


( 6 ) 


- 2(r cose-t)(r-t cose) 2 -> 
(r^+t* - 2rt cose)* J 


(7) 
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With an eye on combining the infinite plane solution with the ring 
solution, we now express the normal and shear components of the stress 
state along the circles r=a and r=b in the plane in terms of the follow’ 
ing Fourier series: 

a lro Ca ' 0 > 'iFtlferf A n (t) Sln n9> 

(8) 

“lrr< a - 9 > * - B n (t) cos ne - 

(9) 

■W»- e > = ^nfcrf c n (t) sinn9 - 

(10) 

W 5 ’ 9 ) * - D " (t) C0s n9 ' 

01) 

where the Fourier coefficients are given by 


„ fJ .s _ 2 f n r 2 sine(2t cose-a) 

A n (t) ‘ ir J C a^+F-'^Ft'Tose 
0 


- W? C0S6? ] S1 " " 9 de - 

(12) 

W { vif ) TT J CT i rr Ca»e)d 0 , 

0 


B n (t) - (- -Lj ) ) ^ j a lrr (a,e)cos ne de; 

0 

(13a, b) 

c n (t) = (~^~f h f a i r0 ( b »e)sin ne de, 

*0 

(14) 

D 0 (t) - (l 2^) ] j o lrr (b,e)de, 

J o 

V*) " ~^f it j a i rr (b> 0 )cos ne de 

^ 0 

(15a ,b) 
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After some manipulations and combining some of the terms, 
grals which appear in (12-15) can be evaluated in closed form, 

i ? a n+1 a n-3 a n +l 

- X { fV - Cn-4)(|-) ]*2n(f) 

o u , n-5 , n-1 

-X&PVU f> -(f) 


B 0 (t> ■ - f * 

B n( t ) = f ( (4-a' i ) 2 ' [‘I’ " f^t) " tf-t' 1 ) 


n+1 


t 2 +a 2 


> n - 3 l,a, n -). l,a 1 l"- 21 - 3 

4 

i 

n+1 


♦ 2(f) fz=r - t"(f) 


n-1 


n-2 | —1 


-(2n-8)(|) + (| n-2 1 -2) (2-)' ]}, n>l, 

i 9 . 4 + n+1 * n+5 + n-1 

V*) ■ F ( TF?FF [ ( F> (f> 

+ ("(f)"" 1 - ( n+4 )(f) n+3 ]1 ’ " ^ 2> 


D 0 (t) 


2t 

w 


0 n (t) . ^ <- ^ [(n*3)(|) 


" +1 n+4/t» n+3 

- inF 


, + I n-2 | +1 2 t n-1 

- £<|n-2|+2)(£) ] + fV<F> 

„. 4 a. n+3 t n+5 , . | n-2 1 +3 

- Xb^FF [( ¥ ) “ ]>. n > 1 , 


b ) The /iXng &ohjutlon. 

In this case assuming a stress function of the form 
^ 2 (r,e) = F(r)e &e 


the inte- 
giving 


(16) 


(17a, b) 


(18) 


(19a ,b) 


(20) 
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the general elasticity solution in polar coordinates may be obtained 
as [9] 


* 2 (r.e) = a 0 + b o log r + c o r2 + d o r2 log r 

+ (a^ + b^ log r + c' Q r 2 + d^r 2 log r)e 

+ (a, r + b,r log r + ^-+ d^ 3 ) 


+ (ajr + b]r log r)e 


+ \ ^V" + V"* + V"" + V*"> cos ne 
n=d 


(21) 


This is basically the solution given by Michell, except that he omitted 
some of the repeated roots. From (21) the stresses for an annulus are 
found to be 


2c l 


_ u o . ~ , , ‘- u i . 9ri sine 

°2rr - r* 2c o + ~F~ + 2d l r) cose 


n-2 


- 1^- (R y sine + R x cose) - z [a n n(n-l)r 
+ b n (n+l)(n-2)r n + c n n(n+l )r"^ n+2 ^ + d n (n-l )(n+2)r" n ]*™ "J , 


2c l 


°2m ~~ F + 2c 0 + + 6d l r >c«e + fe ( V ine + R X C0S9) 

+ z [a n(n-l )r n-2 + b (n+l)(n+2)r n 
n=2 n n 

+ c n n(n+l)r"^ n+2 ^ + d n (n-2) (n-1 )r n ]*™ JJ| , 
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(22a-c) 


_ a 0 J. , 2c 1 _tCOS0 1-V 


a 0 a = dr + (— r - 2d,r) - 7 — (R^cose + R v sine) 
2r8 'r- v r d 1 -sine 4irr y x 


- e [a n n(n-l )r n “ 2 + b n(n+l)r n - c n(n+l)r~^ n+2 ^ 


n=2 


a n ( n n -n-,cos ne 
- d n n(n-l)r J_ s1n n0 » 


where R x and R y are the resultants of the tractions on r=a. In the 
problem under consideration R =0, R =0 and the stresses satisfy the 

X j 

following symmetry conditions: 


a 2rr Cr ’ 0) = a 2rr (r ’” e) ’ a 2ee (r ’ 0) = a ?AB (r ’“ 0) ’ 


2ee' 


a 2re^ r ’ 0 ^ " ' a 2re^ r,_0 ^ 


(23a-c) 


Thus, from (22) and (23) it follows that 


2ci 


n-2 


a 2rr = + 2c Q + (- -p- + 2d 1 r)cose - z [a n n(n-l)r 

+ b n (n+l)(n-2)r n + c n n(n+l )r“ n " 2 + d n (n-l)(n+2)r _11 ]cos ne, 

°2re = (’ 7T- + 2d ] r)sine + s [a n n(n-l )r n “ 2 

+ b n n(n+l)r n - c n n(n+l )r' n “ 2 - d n n(n-l )r _n ]sin ne. 


2ee 


+ 2c + (^- + 6d-,r)cose + e [a n(n-l)r n “ 2 

or I o'* 


+ b n (n+l)(n+2)r n + c n n(n+l )r“ n " 2 + d n (n-2) (n-1 )r“ n ]cos ne . 

(24a-c) 


c) Boundary conditio n& . 

After determining the basic form of the solutions for the infinite 
plate with a dislocation and for the ring, the stresses in the ring 
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having a dislocation may be expressed as the sum of the two solutions 
(see (1), (5-7) and (24)). The combined stress state must then satisfy 
the following boundary conditions (Figure la): 


a^ rQ (a,Q) 

°lrr^*® ) 

tf lr6 ( b »0) 

a le8 (r, 0 ) 


+ <J 2 rr (a,e) 

+ *2r0 (a ’ 0) 
+ a 2rr^’ 9 ^ 

+ * 2r0 (b,0) 
+ cr 200 (r,O) 


= 0, 0 <_ 0 <_ tr , 

= 0, 0 <_ 0 <_ u , 

= 0, 0 £ 9 <_ ir , 

= 0, 0 <_ 9 <_ TT 

= g(r), c < r < d 


(25 a-d) 
(26) 


where a-| rr and a^ r0 are given by (8-15) and a 2rr and a 2r0 are given by 
(24), and g(r) is the crack surface traction. Note that each equation 
in (25) is a sine or cosine series in which the coefficients of sin ne 
and cos ne, (n=Q,l,...) must vanish. Thus, defining 


a n = “ 


uf 

■jr(ic+l ) a r 


uf 

ir(ic+l ) 




c n 


uf 

TT (k+1 ) 


uf * 
tt(k+ 1 ) n 


(n— 0 ,1,...) , 


we obtain 


(27 a-d) 


3 3 

# + 2y o = ' B o ’ W * ^ 0 * -°o ; 


2yi 2yi 

- -pL + 25l a - - Bl , + 25-jb = -D-j ; 


2Y1 


+ 26 


i a ' A i • - + 


26^ = C 1 ; 


(28a, b) 
(29a, b) 
(30a, b) 
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<* n n(n-l)a n " 2 + 3 n (n+l )(n-2)a n + v n n(n+l)a n ’ 2 
+ S n (n-1)(n+2)a~ n = B n , (n=2,3,...) , 
a n n(n-l)a n-2 + 8 n n(n+1)a n - y n n(n+l)a" n ” 2 

- <5^ n(n-l )a - , (n = 2,3,...). , 

a n n(n-l)b n “ 2 + s n (n+l)(n-2)b n + Y n n(n+l)b' n " 2 
+ 5 n (n-l)(n+2)b" n * D n , (n=2,3,...), 
a n n(n-l)b n ~ 2 + s n n(n+l )b n - y n n(n+l)b" n “ 2 

- 5 n n(n-l)b" n = C n , (n=2,3,...) , (3la-d) 


where A n , B n , C n , and D n , (n=0,l ,2, . . . ) are known functions of t and 
are given by (16-19). Consequently, the coefficients a n , b n , c n , and 
d n calculated from (27-31) will also be functions of t. From (24) it 
may be observed that the coefficients a Q , a-| , b^ and d Q do not appear 
in the expressions of ^ij’ Therefore, if the series in (24) are 
truncated at the nth term, there would be 4n unknown coefficients to be 
determined. Cn the other hand, there are 4n+2 equations in the corre- 
sponding algebraic system given by (28-31). Two of these equations must, 
therefore, be redundant. This is, in fact, suggested by the pairs of 
equations (29a, 30a) and (29b, 30b). These equations indicate that if 

A-j + B-, = 0 , C-, + D-| = 0 , (32a, b) 

then (29) and (30) are identical and by, for example, ignoring (30), 
the unknown coefficients may be determined uniquely from (28), (29), 
and (31). Referring to (8) and (9), it may be seen that 
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2v 

A-j + b t = f [a lrQ (a,e}sin 0 - <j, (a,0)cose]ade. 

(33) 

Consider now an infinite plane containing an edge dislocation at x=t 
(see eqs. (2,3)). If we introduce a circular hole r=a, a<t into the 
plane (Fig. 1), then certain tractions a^ rr and a-| r0 must be applied 
along the boundary of the hole in order to maintain -the equilibrium of 
the plane. It is clear that the integral in (33) is nothing but the 
x-component of the resultant of these tractions. Since the external 
force system due to the dislocation or the crack (which results from 
(2) and (3)) is self-equilibrating, this resultant must be zero, proving 
the validity of (32a). Equation (32b) follows from similar arguments 
for a disk r<b containing a dislocation at x=t<b. It is because of the 
relations (32) that in (16-19) the expressions for A-j and C-| are not 
included. 

After expressing the coefficients a n d n in terms of f, the 

remaining boundary condition (26) may be used to determine f. If we 
now assume that the crack is formed by distributing the dislocations 
along the line 0=0, c<r<d with f as a function of t, o Qa (r,e) in the 
ring may be evaluated by integrating the Green's function obtained from 
(1), (7), (24c) and (27-32). Specifically, for 0=0 from (26) we find 


’ 


lit! 

t-r 


dt + 


d 

‘ 

■ 

c 


k(r,t)f(t)dt 


tt(k+1 ) 
2y 


g(r). 


c < r < d, 


(34) 


where 

k(r,t) = “ \ pr + 2 y 0 + ~^r + 66^r + r [n(n-l )<* n r n " 2 
+ (n+l)(n+2)B n r n + n(n+l )y n r" n ”^ 

+ (n-2)(n-l)6 n r’ n ], (35) 
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and o n> e n> Y n > and 6 n (n=0,l,...) are known functions of t which are 
determined from (28-32). If the crack is an embedded crack (i.e., if 
a<c<d<b) then from the definition of f(t) given by (3) it follows that 



f(t)dt = 0 


C36) 


3. Stress Intensity Factors and Examples 

If the crack is an embedded crack, the index of the integral equa 
tion (34) is +1 and the solution is of the following form 

f(r) = L(r-ci&r)Ji ’ (c<r<d) (37) 

where F(r) is a bounded function. After evaluating F(r), the stress 
intensity factors at the crack tips r=c and r=d may be defined and 
obtained from 


k(c) = lim /2(c-r) a 0Q (r,0) 
r>c 


= i im v^[r^cT f (r) 
r-x: 


2u F(c) 

1+K /(d-c)/2 ’ 


k(d) = lim /2(r-d) a (r,0) 
r-xJ 


2p_ 

1+K 


lim /2(d-r) f (r) = 
r-xi 


2u F(d) 
1+K /(d-c)/2 


(38) 


(39) 


The singular integral equation (34) may be solved numerically by 
normalizing the interval (c,d) to (-1,1) and by using the Gaussian inte- 
gration formula [10], Since the Fredholm kernel k(r,t) is given in a 
simple form, one could obtain the results to any desired degree of 
accuracy with a relatively modest computational effort. 

In the case of an edge crack (i.e., for a=c<d<b or a<c<d=b) the 
single-valuedness condition (36) is not valid and also is not required 
as an additional condition for a unique solution. 
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For certain ring-crack geometries the convergence of the infinite 
series giving the Fredholm kernel (35) was rather slow. The results 
given in this paper were obtained by truncating the series at the Nth 
term, solving the integral equations for N=100, 120, and 140 for each 
ring-crack geometry and loading condition considered, calculating the 
stress intensity factors k as a function of N (from C38) and (39)), 
and by using a curve-fitting procedure of the following form 

k(Nj) = A + JL , (N 1 = 100, N 2 = 120, N 3 = 140) . (40) 

The value of A obtained from (40 ) is assumed to be the calculated 
stress intensity factor. In the examples, for the values considered, 
a is generally greater than unity, implying that the slope of k vs. 1/N 
curve at (1/N) = 0 is zero. Also, since the values k(N.) (i=l,2,3) 
differ from each other only in the third or fourth significant digits 
in the slowest converging cases, the extrapolated values ought to be 
quite reliable. 

Three different loading conditions are used in the examples. Note 
that g(r) which appear as the input function in the integral equation 
(34) or as the traction in the boundary condition (26) is the value of 
a Q0 in the perturbation problem. For the loading conditions considered, 
g(r) is given by the following expressions: 

a\ Uni^oAm osumk Au/ifiace. psizAAuAZ 

g(r) = -ct q , c < r < d. (41) 

b) Inte/inatty pH.2A6usu.zzd cyUnd&H. 

g(r) = " g 2 r a 2 ~ ,c<r<d, (42) 

where p Q is the internal pressure (i.e., a rr (a,e) = -P 0 ). Here, it 
should be observed that in the case of the internal edge crack (i.e.. 
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for c = a < d < b), in addition to g(r) given by (42), the crack surfaces 
may be subjected to uniform (fluid) pressure g(r) = -p Q . If that is the 
case, then the stress intensity factor separately obtained from gCr) = 

-p Q should be added to that given by (42). In the examples, results 
due to (41) and (42) are listed separately. 

c) Rotating diik 

g(r) - - ^ p [a* + b* + (^) 2 - ^ r*] (43) 

where p is the mass density and u is the angular velocity. The expression 
given in (43) is for the generalized plane stress and is valid for a 
rotating "disk". For a long rotating hollow cylinder, i.e., for the 
plane strain case, we have 

9(r) = " i T^T p [ b2 + a2 + ( 7 ^) " 3127 r 2 ]. (44) 

For various crack geometries the calculated stress intensity fac- 
tors are shown in Tables 1-8. To give some idea about the trends, some 
limited results are also given in Figures 2-5. Table 1 gives the results 
for an embedded crack symmetrically located with respect to the boun- 
daries of the ring or the hollow cylinder and subjected to a uniform 
crack surface pressure (see (41)). Here the crack length is fixed at 
half the wall thickness b-a and the ratio of a/b is varied. The results 
are also displayed in Figure 2. Intuitively, it may be argued that 
the "constraint" at the inner crack tip r=c is greater and at the outer 
crack tip r=d is less than that of a centrally cracked infinite strip. 
Therefore, one may expect that k(c) should be less and k(d) should be 
greater than the corresponding stress intensity factor in the strip. 

This trend may also be observed in Figure 2 where the ring results seem 
to approach the strip value as a-*» for constant wall thickness b-a. 

Table 2 shows the effect of the crack length on the stress inten- 
sity factors in a ring or a hollow cylinder with a fixed radius ratio 
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Table 1. 


Stress intensity factors k(c) and k(d) for an 
embedded crack in a ring or a hollow cylinder sub- 
jected to uniform crack surface pressure a QO (r,0)= 
d-c _ - - 09 


V 


Fa 


= 0.5, 


c-a _ 
F?" 


1 (Figure la) 


a 

Fa 

k(c) 

kfd) 

V(d-cl/2 

<V(d-c)/2 

0.05 

1.1477 

1.2046 

0.10 

1.1498 

1.2030 

0.25 

1.1580 

1.2018 

0.50 

1.1664 

1.2007 

1.0 

1.1736 

1.1980 

2.0 

1.1788 

1.1943 

3.0 

1.1810 

1.1923 

4.0 

1.1822 

1.1911 

CO 

1.1867 

1.1867 


Table 2. Stress intensity factors k(c) and k(d) for an 
embedded crack in a ring or a hollow cylinder 
subjected to uniform crack surface pressure 
d 0e (r,O) = -a Q (columns 2 and 3) or internal 
pressure a rr (a,e) = -p 0 (columns 4 and 5). 

■| = p = 1 (figure la) (The last column 

refers to the centrally cracked infinite strip) 


d-c 

Fa 

Me) a 

kfd) 

k(c) 

k(d) 

k 

cr 0 /[t‘C)/Z 

V(d-c)/2 

P</(d-c)/2 

P 0 /(d-c)/2 

<V(d-c)/2 

+0.0 

1.0 

1.0 

0.9259 

0.9259 

1.0 

0.1 

1.0060 

1.0062 

0.9522 

0.9129 

1.0060 

0.2 

1.0241 

1.0256 

0.9920 

0.9139 

1.0246 

0.3 

1.0556 

1.0608 

1.0468 

0.9311 

1.0577 

0.4 

1.1034 

1.1158 

1.1201 

0.9682 

1.1094 

0.5 

1.1736 

1.1980 

1.2186 

1.0318 

1.1867 

0.6 

1.2785 

1.3213 

1.3557 

1.1341 

1.3033 

0.7 

1.4445 

1.5146 

1.5606 

1.3010 

1.4884 

0.8 

1.7407 

1.8536 

1.9084 

1.6002 

1.8169 

+1.0 

00 

00 

00 

00 

00 
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Table 3. Stress intensity factors k(c) and k(d) for an 
embedded crack in a ring or a hollow cylinder 
subjected to uniform crack surface pressure 
o 0e (r,O) = -cr 0 (columns 2 and 3) or internal 
pressure a rr ( a,e) = -p 0 (columns 4 and 5). 

0.5, |= \ (Figure la) 


c-a 

k(c) 

k(d) 

k(c) 

k(d) 

o 0 /(d-c)/2 

V(d-c)/2 

Pc/(d-c)/2 

P</ld-cJ/7 

0.15 

1.2414 

1.1755 

1.4292 

1.1004 

0.20 

1.1929 

1.1777 

1.3033 

1.0554 

0.25 

1.1736 

1.1980 

1.2186 

1.0318 

0.30 

1.1744 

1.2391 

1.1608 

1 .0290 

0.35 

1.1936 

1.3118 

1.1249 

1.0533 


a/b = 1/2. The crack is again synrnetrically located with respect to 
the boundaries (i.e., c-a = b-d). The table shows the results for uni- 
form crack surface pressure o q (columns 2 and 3), for internal (wall) 
pressure p Q (columns 4 and 5), and, for the purpose of comparison, for 
a centrally cracked infinite strip (column 6). Again to show the gen- 
eral trend of the results, the stress intensity factors for a cylinder 
under internal pressure p 0 are also displayed in Figure 3. In this 
case, as the crack length d-c goes to zero, the stress intensity fac- 
tors k(c) and k(d) approach a 0Q ((^£), 0)/(d-e)/2 which is the corre- 
sponding infinite plane result (for a/b=l/2, o 00 ((^-, 0) = (25/27) p Q = 
0.9259 p 0 ). In this problem the initial decrease in k(d) with increas- 
ing crack length is due to the decrease in cr 06 (r,O) in the neighborhood 
of r=d. 

The effect of the relative position of the crack in the cylinder 
wall is shown in Table 3. Here it is assumed that the crack length and 
radius ratio of the cylinder are constant ((d-c) = (b-a)/2, a = b/2) 
and the radial position of the crack (c-a)/(b-a) is varied. 

The results for an internal edge crack c = a<d<b are given in 
Tables 4-6. Table 4 shows the stress intensity factor for a uniform 
crack surface pressure a Q (see (41)). The strip results given in the 
last column of the table are the limit of the cylinder results for 
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Table 4. The normalized stress intensity factor k(d)/(<r 0 /d-a) 
in a hollow cylinder or a ring containing an internal 
edge crack (c=a<d<b) and subjected to uniform crack 
surface pressure cr aQ (r,0) = -a . 

u u O 


v a 



d-X 5 ^ 
FT \ 

1/3 

1/2 

->0 

1.1216 

1.1216 

0.1 

1.153 

1.155 

0.2 

1.218 

1.229 

0.3 

1.295 

1.310 

0.4 

1.373 

1.402 

0.5 

1.465 

1.508 

0.6 

1.578 

1.635 

0.7 

1.730 

1.796 

-*■1 . 0 

00 

00 


1 

2 

3 

00 

(strip) 

1.1216 

1.1216 

1.1216 

1.1216 

1.157 

1.159 

1.167 

1.1893 

1.247 

1.277 

1.299 

1.3674 

1.366 

1.449 

1.493 

1.6601 

1.503 

1.655 

1.747 

2.1119 

1.658 

1.901 

2.066 

2.8258 

1.830 

2.177 

2.441 

4.035 

2.030 

2.475 

2.851 

6.361 

CO 

CO 

oo 

00 


Table 5. 


The normalized stress intensity factor k(d)/[a 08 (a,O)/d-a] 
for an internal edge crack (c=a<d<b) in a holl&w cylin- 


der under internal pressure, a rr (a,e) = -Pg. (The effect 
of pressure p on the crack surfaces are not included) 
a 0e (a,O) = P 0 (b 2 +a 2 )/(b 2 -a 2 ). 


a 

. X E-a 

d-a \ 
b-a 

1/3 

1/2 

1 

2 

3 

+0 

1.1216 

1.1216 

1.1216 

1.1216 

1.1216 

0.1 

0.878 

0.957 

1.058 

1.113 

1.137 

0.2 

0.756 

0.878 

1.057 

1.189 

1.239 

0.3 

0.685 

0.834 

1.089 

1.305 

1.397 

0.4 

0.643 

0.813 

1.138 

1.461 

1.613 

0.5 

0.620 

0.808 

1.199 

1.635 

1.883 

0.6 

0.611 

0.816 

1.273 

1.834 

2.189 

0.7 

0.616 

0.841 

1.359 

2.046 

2.527 

+1 

00 

00 

00 

00 

oo 
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Table 6. 


The normalized stress intensity factor k(d)/[agg(a,0)/d-a] 
for an internal edge crack (c=a<d<b) in a rotating disk 


(v=0.3) a flfl Ca,0) = 3+v 

4 


p or 


( 5 2 


+ 1 -v a' 
3+v 




a 

IPa \ 

1/3 

1/2 

1 

2 

3 

->0.0 

1.1216 

1.1216 

1.1216 

1.1216 

1.1216 

0.1 

1.008 

1.040 

1.089 

1.120 

1.139 

0.2 

0.962 

1.018 

1.114 

1.202 

1.240 

0.3 

0.944 

1.018 

1.168 

1.324 

1.401 

0.4 

0.943 

1.032 

1.235 

1.477 

1.617 

0.5 

0.956 

1.058 

1.315 

1.664 

1.879 

0.6 

0.981 

1.096 

1.405 

1.867 

2.196 

0.7 

1.025 

1.153 

1.508 

2.082 

2.526 

->1.0 

CO 

CO 

00 

oo 

00 


a-*» and for constant wall thickness b-a. The results of Table 4 are 
partly displayed in Figure 4 to show the general trend. One may again 
note that as a/(b-a) increases, the geometric constraints decrease caus- 
ing an increase in the stress intensity factors. As the crack length 
d-a tends to zero, the stress intensity factor ratio approaches that 
of a half plane with an edge crack namely, 1.1216. 

Table 5 and Figure 5 show the stress intensity factor k(d) for 
a hollow cylinder with an internal edge crack and subjected to internal 
pressure a rr (a,e) = -p Q . Note that k(d) is normalized with respect to 
the hoop stress at the inner boundary of the cylinder, a 00 (a,e) = 
P 0 (b 2 +a 2 )/(b 2 -a 2 ). Figure 5 also shows the variation of the hoop stress 
dee (r,e) through the cylinder wall, again normalized with respect to 
a Q0 (a,e). Thus, the initial decrease in k(d) for increasing crack length 
d-a in cylinders having small values of a/(b-a) may easily be explained 
by the relatively sharp decrease in a 00 (r,e) (which is the external load 
applied to the crack surface). Again, as d-a->0, the stress intensity 
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factor ratio tends to the half-plane value 1.1216. 

The rotating disk or cylinder results are given in Table 6. Here, 
too, the stress intensity factor k(d) is normalized with respect to 
the hoop stress a 0Q at the inner boundary r=a, the expression of which 
is given in the table (see (43)). Note that the expressions of the 
hoop stress c Qe in the rotating disk and in the long hollow cylinder 
are different and are dependent on the Poisson's ratio v. Therefore, 
the problems for the disk and the cylinder must be solved separately for 
a specified value of v. The results given in Tables 6 and 8 are for 
the rotating disk which are obtained by using the plane stress solution 
with v=0.3. 

The results for a hollow cylinder or a disk with an external edge 
crack are given in Tables 7 and 8. Table 7 shows the results for an 
internally pressurized hollow cylinder or disk. The stress intensity 
factor for a rotating disk obtained again from the plane stress solution 
with v=0.3 is given in Table 8. 


Table 7. The normalized stress intensity factor k(c)/[CT 00 (b,O)> / b-c] 
for an external edge crack (a<c<d=b) in a hollow cylinder 
or a dtsk under internal pressure a rr (a,e) = -p 0 ; cr 00 (b,O) 
2a 2 p 0 /(b 2 -a 2 ). 


a 

. b-a 
b-cy 
b-a 

1/3 

1/2 

1 

2 

3 

■*0 

1.1216 

1.1216 

1.1216 

1.1216 

1.1216 

0.1 

1.250 

1.245 

1.219 

1.203 

1.208 

0.2 

1.413 

1.391 

1.368 

1.355 

1.352 

0.3 

1.617 

1.587 

1.557 

1.560 

1.574 

0.4 

1.874 

1.829 

1.797 

1.832 

1.873 

0.5 

2.207 

2.137 

2.095 

2.170 

2.259 

0.6 

2.655 

2.540 

2.462 

2.586 

2.740 

0.7 

3.301 

3.103 

2.941 

3.085 

3.320 

+ 1.0 

-*» 



-*» 

-*» 
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Table 8. The normalized stress intensity factor k(c)/(a 0 n(5,O)/b^c) 
in a rotating disk having an external edge cracx (a<c<d=b) 
a ftfi (b,0) = 3+v p ui 2 (a 2 + 1-v b 2 ), v = 0.3. 

~T~ 3+v 


\ a 

\ b-a 
b-c\ 
b-a 

1/3 

1/2 

1 

2 

3 

+0 

1.1216 

1.1216 

1.1216 

1.1216 

1.1216 

0.1 

1.312 

1.288 

1.245 

1.217 

1.214 

0.2 

1.530 

1.486 

1.419 

1.377 

1.367 

0.3 

1.774 

1.714 

1.632 

1.608 

1.601 

0.4 

2.053 

1.982 

1.893 

1.894 

1.911 

0.5 

2.380 

2.299 

2.210 

2.242 

2.309 

0.6 

2.777 

2.691 

2.601 

2.682 

2.811 

0.7 

3.301 

3.210 

3.093 

3.199 

3.406 

+1.0 

-x» 


-*co 

«*x> 
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Figure 








(d-c) / (b-a) 


Figure 3. 


Stress intensity factors for an embedded crack in a hollow 
cylinder or a disk under internal pressure a rr (a,0) = -p Q . 






Figure 4. The stress intensity factor for an internal edge crack in 
a hollow cylinder or a disk loaded by a uniform crack sur- 
face pressure a 00 (r,O) = -a Q , £=d-a. 



Figure 5. The stress intensity factor ratio k' = k(d)/[a 0 ^(a,O)/d-a] 
and the normalized hoop stress a 0 g(r)/o 00 (a,O) in a hollow 
cylinder or a disk containing an internal edge crack and 
subjected to the internal pressure a rr (a,8) = -p Q ; cr 00 (a,O) = 
p 0 (b 2 +a 2 )/(b 2 -a 2 ) (The effect of the pressure p 0 on tne 
crack surfaces are not included). 
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